Quantum Phases of Bose-Hubbard Model in Optical Superlattices 
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In this paper, we analyze the quantum phases of multiple component Bose-Hubbard model in 
optical superlattices, using a mean-field method, the decoupling approximation. We find that the 
phase diagrams exhibit complected patterns and regions with various Charge Density Wave (CDW) 
for both one- and two- component cases. We also analyze the effective spin dynamics for the two- 
component case in strong-coupling region at unit filling, and show the possible existence of a Spin 
Density Wave (SDW) order. 
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I. INTRODUCTION 

The observation of Mott insulator - superfluid tran- 
sition of ultracold bosons loaded in optical lattices [l| 
has triggered huge amount of interest in quantum 
simulation [2]. A lot of efforts have been made to in- 
vestigate Bose-Hubbard Hamiltonian both sing lei 
and multiple components Q, in cold atom systems 
(in double- wells Q, superlattices 0, Q)^ using vari- 
ous techniques ( proj ection wave- funct ion 0, decoupling 
approximation pi Hoj. field theory [111], dynamical mean- 
field [l2l]. etc.), aiming to achieve a comprehensible un- 
derstanding of this many-body model 13]. 

Among these, multi-component Bose-Hubbard in op- 
tical superlattices attracts special attention due to its 
rich quantum phases [3, [3 and promising potential in 
emulating effective spin-spin interactions @, Ha]- Many 
works on this topic have revealed the existence of Mott- 
insulator with integer and fractional fillings, which cor- 
respond to various Charge Density Wave (CDW) and 
complected dynamical behaviors [lE d, E, HI H3 • 
People also reported topological phases and exotic lo- 
calizations in disordered superlattices [2 U l24j. Recently, 
to obtain a complete and thorough description of this 
model, several attempts have been made pa. |26|. 

In this paper, we reinvestigate this problem by us- 
ing a mean-field approach, the decoupling approxima- 
tion. Working in the Mott insulating phase as an unper- 
turbed ground state, and treating the hopping as pertur- 
bations, this method (though at mean-field level) can still 
give comparatively accurate results for phase transition, 
comparing with numerical simulations [4|. Besides, since 
this method decouples the original Hamiltonian into a set 
of single-site Hamiltonian, it exhibits simple energy ex- 
pressions for different filling situations in the Mott state. 
Therefore, we can analyze and write down all the possi- 
ble filling configurations in a clear and evident manner. 
Then, We depict the corresponding phase boundaries for 



each filling condition, and gain a systematic description 
of possible quantum phases of this model. We further 
consider the effective spin dynamics at strong coupling 
limit and unit filling, and find both a ferromagnetic and 
an anti-ferromagnetic spin- wave excitation as the poten- 
tial barrier between neighboring sites varies. 

The paper is organized as follow: We first introduce 
the model in Sec II, then deals with the single-component 
Bose-Hubbard model in Sec. Ill, which acts as a prece- 
dent of the two-component case that is discussed in Sec. 
IV. Both sections contain the mean-field calculations, the 
analysis of possible filling patterns and phase diagrams 
with corresponding illustrations. The spin dynamics is 
analyzed in the last subsection in Sec. IV. We give the 
conclusion in Sec. V. 



II. MODELS 

Our starting point is a two-component Bose-Hubbard 
model in a double-periodic superlattices, 

H = - {Uja^aja + h.c.) + — ^ h lrT (h ia - 1) 



(1) 
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Here (ij) denotes the nearest-neighbor counting, is the 
hopping amplitude, U is the inter-species repulsion, V is 
the intra-species repulsion, [i a is the chemical potential 
that restricts the particle number, Aj is the energy bias 
on a given site, hi„ — al a ai a is the number operator for 
bosons with a =f , J. representing the two internal states 
of the trapped Bose Einstein Condensate (BEC). 

For single-component Bose-Hubbard, the (pseudo)spin 
gets polarized so that a takes a certain value; thus there 
is no intra-species interaction V . 

Generally speaking, the period of the superlattice can 
be set to be an arbitrary integer /, as long as we require 
physical quantities are periodic functions, ipi = For 



2 



one supercell^ 




FIG. 1: (color online) An illustration of a period-2 super- 
lattice. The potential barrier is A. 1 (2) is denoted as the 
deep (shallow) site in a supercell. The red and blue arrows 
represent different internal states of the trapped BEC. 

I = 2, we can simply set the potential bias to be 

f o,i G odd (2) 

[ A,i€ even ' v ; 

as shown in Fig. fl}. 



III. ONE-COMPONENT CASE 

A. Basic Formulism 

First, let us consider the one-component case for clar- 
ity. The Hamiltonian reduces to 

H = - ^{tijal&j + h.c.) + ^J2ni(hi~ 1) 
(«'} » 
-^(M-Ai)ni. (3) 

i 

To determine the phase boundary of Mott-insulator (MI) 
- superfluid (SF) transition, we apply an extended decou- 
pling approximation which was developed in Ref. Q . In 
this mean-field scheme, we write the two bosonic opera- 
tors a\a,j as a\dj « (o-l)&j + — (&J )(%)• Then by 
introducing a site-dependent (local and small) SF order 
parameter, ipi = (a|) = (<Zj), the hopping term can be 
written as, 

_ 53(*y fi *Sj + ft-c) » -ty^[^i(aj + a + aj +5 ) - ipiip i+s 

(ij) i,8 

+ip l+ s(a l + a}) - ipi+sipi] 
= -2t y^(a| + eg - ipi) ^ 

i 5 

= -2*t^(at+Oi-^i)V>i+i, (4) 

where we assume that in the strong-coupling regime 
(U, V tij, tji) tij m tji = t and denote ipi+s — 4>i-s 

5 is the position vector for the nearest neighbors, 



z is the coordination number. We then decouple the ini- 
tial Hamiltonian Eq. © into N (nearly) independent 

parts 0], H = Yli=i Hit where 

Hi = (hi - 1) - (/x - Aj) h t 

-2zt(a\ + &i - ipi)ip i+ i. (5) 

As a result, we can concentrate in one super-cell as shown 
in Fig. {T]) with reduced Hamiltonian 

2 

H sc = {2zty 1 Y J Hi 
l=i 

2 U 
/=i 

-(a| + ai - in)i]ji + i] 
= i> 2 V x +^V 2 , (6) 

where the dimensionless quantities are U = U /2zt, p, = 
n/2zt and A/ = Ai/2zt, the unperturbed Hamiltonian 

m is 

= E["2 ( "< " 1} " ^ ] + 2M2 (7) 

with reduced chemical potential /2; = /2 — A; , while the 
hopping-induced perturbations are 

V, = -(hj + &i). (8) 

In deriving above equations, we have used the periodic 
condition of the superlattice, V'/+2 = 4>h I — lj 2. 

The unperturbed energy is the energy of states that 
have definite particle number (Mott state), namely, 

= ^. 92} = min{^)}„ =0il , 2 ,..., (9) 

where g = gi + .92 is the average particle number in one 
supercell, gi is the particle number on a certain site I in 
the supercell. Up to the second order, this implies 

f <*,f (10) 

thus there will be a constraint on the chemical potential 
ft. (To be discussed later.) 

In the Mott state and near the MI-SF transition, the 
on-site particle number is still well-defined, then the un- 
perturbed energy (written in a dimensionless form) be- 
comes 

2 fj 

4L 2 } = E^-5/(2i-l)-MiS/] + 2V^2 

= \ (gi + 92 -g)- fig + A.92 + 2^1^11) 

where we have used Eq. 
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Meanwhile, the perturbation term V results in an en- 
ergy correction, we can calculate this second-order per- 

(2) r — I 

turbation energy Eg in a standard manner [J] , 



E { 2) = Ef ] , 
s {91-92} 



v4 



E 



\(ni;n 2 \V 1 \g 1 ;g 2 )\ 2 



E 



(0) 



E 



(0) 



0? 



ni+"2#9l+92 ^{51:52} "{ni;n 2 } 

|(ni;n 2 |y 2 |5i;52)| 2 



E 



E 



(o) 



™i+n2#9i+92 {91:92} 



(0) 

{m;n.2} 



(12) 



where the summation over {774 ; n 2 } has only the following 
terms, {n\\ ra 2 } = {51 ± 1; 52} , {51; 32 ± 1}, due to the 
simple form of Vi . After straightforward calculations, we 
find 



E 



(2) 

{91:92} 



9i 



1 



9i 



-Ugt+p U( gi -1) 
92 + 1 



+ 



/' 

.92 



-Ug 2 +p-A U(g 2 -l)-p, + A 
-i>\ (p + U) 



[U( 9l -l)-p] {-Ugx+p) 

-tf jfi ~ A + g) 

[^(fl2 — 1) — M + A] (-Ug 2 +p- 

Therefore, the total energy of a super-cell is E = E, 



(13) 



(0) 



t 2 \ 

Eg + ■ ■ ■ can be expanded into a power series of the SF 
order parameter ^ and ip 2 (Landau expansion), 

E (ipi,ip2) =a a + a 2 i\)\ + b 2 ipl + c 2 ipi^ 2 + O (ipf, ip 2 ) , 

(14) 

where we presume all coefficients of fourth-order and 
above are positive, stabilizing the system. The coeffi- 
cients are 



«o 



«2 



\ (si + .92 ~ .9) ~ M.9 + A52, 



(p-A + U) 



[U (<7 2 - 1) - p + A] (-Ug 2 + p-A 

-(p + u) 



(15a) 
,(15b) 



[U( gi -l)-p] (-U91+H) 



c 2 = 2. (15c) 



Clearly, E (ipi — 0, ip 2 = 0) is a local extremum, and 
represents the MI state. As the interaction U, the 
biased potential A and the particle number (chemi- 
cal potential p) changes to a critical point, the ex- 
tremum E (0, 0) becomes an instability point, thus 
a phase transition to SF state occurs. To express 
this idea explicitly, we should compare <9l E\ (o,o) > 
d\ £71(0,0) an d ^jfl^jJz/lfOjQ)- The critical condition is 

^^(o.o^^o.o) = i^id^E] (0)0) ) 2 , i.e., 4a 2 6 2 = c?,, 
or equivalently, 



[C7 (. 9l - 1) - A] (-CTpi + M) (92 - 1) - + A] (-C7 52 + /I - A) ■ 
Simplify it, we obtain the phase boundary, 

- = 1 / [( gl — 1) — p] (- 9l + jjj [( g2 - 1) - A + A](- g2 + A - A) 
c 2 V (A + i)(A-A + i) 



where t c = zt/f/, p = p/U, and A = A/U. 

Before we come to the phase diagram, let us make 
some remarks on our approximation method. As pointed 
out in Ref. [l(|, the decoupling scheme works well in 
the vicinity of Mott state. The hopping effect induces a 
weak and local superfluid ipi, serving as a perturbation 
over the Mott state. Consequently, although this tech- 
nique can determine the phase boundary conveniently, 
it cannot extrapolate correct physics in deep superfluid 
phase. Besides, from Eqs. (JH1 El an d the introduc- 
tion of many dimcnsionless quantities we know that the 
MI-SF transition is universal, i.e., it occurs in similar 
manners in different dimensions; and this approximation 
gets better accuracy in higher dimensions. 



B. Phase Diagrams 

To depict the phase diagram, we can plot the phase 
boundary surface, Eq. (| 1 T[) in the p-A-t coordinate sys- 
tem. Before that, we need to take a closer look at the 
unperturbed groundstate to find appropriate constraints 
on the chemical potential p. 

As we mentioned, assuming the average filling in a 
supercell is g = g\ + g 2 = {gi;g 2 }, a natural question 
arises: what is the most energetically favorable filling 
factor {gi;g 2 } for a given gl A straightforward solu- 
tion to this question is to compare two arbitrary fill- 
ing configurations {k\; k 2 } and {mi; m 2 } with k\ + k 2 — 
mi + m 2 — g; namely, to compare their unperturbed en- 
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ergy, AE = - figU^. By using Eq. flTJ, this 

energy difference is 

A£ = (fc 2 + m 2 - g) (fc 2 - m 2 ) + (fc 2 - m 2 ) A (18) 

with AE = AE/U, /c 2 ,m 2 = 0,1,..., g, respectively. 
Therefore, for a given potential bias A = A/U and filling 
number g, we can determine the groundstate configura- 
tion {.9i;.g 2 } which has the largest differences AE with 
respect to all other configurations. 

More intuitively, we can simply draw the energy lev- 



els in the two sites of a supercell with a given reduced 
potential bias A, and then fill atoms from the lowest 
level to higher levels one by one, up to g. This filling se- 
quence naturally costs the least energy, thus the resulting 
configuration is just the groundstate filling {gi;<? 2 } (see 
Appendix) . 

As a result, for a given groundstate {<7i;<? 2 }, from Eq. 
(fTOf we obtain the constraint on chemical potential jl. 
For the most imbalanced case {<?; 0} , p, € (g — 1, g). For 
other cases, the calculations are simple but tedious. The 
results are listed in the following table, (p S even.) 



A 


g G even 


g G odd 


(i>- hp) 


/ ue(i[<7+(p-2)],i(<7-p) + A) 


Me(i[<7-(p + l)]+A,i[ 5 +(p-l)]) 


(p,P+l) 


Ae(^-(p + 2)] + A,i(9+p)) 


Me(i[<7+(p-l)],i[<7-(p+l)] + A) 



0.4 s 



0.3 




FIG. 2: (color online) The phase diagram for one-component 
case. Regions under the lopes are Mott insulating phase (MI) , 
regions above those lopes are superfluid phase. 



With these restrictions in hand, we now plot the phase 
diagram for one-component Bose-Hubbard model in su- 
perlattices, as shown in Figs. @ and ©■ It can be seen 
that in strong-coupling regions (t <C U, under the lobes) 
with given potential bias A, the system is in Mott phase. 
As U decreases, phase transitions occur, the system goes 
into superfluid phase that is above the lobes. On the 
other hand, when A changes, the phase boundaries and 
allowed Mott states change correspondingly. 

For clarity, we denote the groundstate filling config- 
uration at different parameter regions in Fig. ([3]). In 
regions to the left of the black solid line, parameters are 




1 2 3 4 5 



fjAJ 

FIG. 3: (color online) A contour diagram for one-component 
case. Regions with different numbers correspond to different 
unperturbed ground state filling patterns. The left (right) 
number represents the particle number in deep (shallow) sit 
in a supercell [see Fig. Q]. This filling configuration reflects 
a CDW order in the Mott state. Regions to the left of the 
black solid line favor an imbalanced filling pattern {g; 0}. The 
chemical potential /i (particle number g) is fixed along the 
vertical (oblique) black dashed lines. 



chosen such that in Mott phase particles in a supercell 
must reside in the deep site (site 1), thus forms crystalline 
structure that has a doubled period of the original opti- 
cal lattice. To illustrate the effect of increasing potential 
bias, we draw the phase boundaries along the vertical 
black dashed lines with fixed chemical potential /i, and 
along the oblique black dashed lines with fixed particle 
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FIG. 4: The phase boundaries along the dashed lines in Fig. 
0. In (a) and (b) chemical potential /i is fixed, in (c) and 
(d) particle number g is fixed. Notations are defined in Figs. 
(|2| and ((3|. SF stands for superfluid phase. 



numbers g in Figs. 



IV. TWO-COMPONENT CASE 



The effective Hamiltonian in one supercell H sc can still be 
written into two parts, H sc — (2zt) 1 J2i=i Hi = + 
V, where the unperturbed part is 



U 

H^) = X! hr"J» (" icr _ !) _ W^^ct + 1plaiH+l,a] 



1=1, 



1=1 



(21) 



with fiiu = ficj — Ai and = V/2zt; while the perturba- 
tion is 



v=-J2 + = v^i+i,' ( 22 ) 



i=i,, 



/=!,. 



with = -(aJ CT + a iff ). 

Similar to the one-component case, the unperturbed 
energy for a supercell with g atoms is = 



i?[ '' , = minion"'' In— i 2 i which means 

E g ^ < Eg^. Thus, the zero-order energy is 



,(o)- 



A. Basic formulism 

The Hamiltonian is just Eq. (fTJ). The hopping term 
can be similarly decoupled into a summation over single- 
site parts, a\ a a jry = (atjojv + a\ a {a jtJ ) - (at 7 )(a i<7 ), as- 
suming the SF order parameter is also spin-dependent, 
ijjia- = {a\ a ) = (ai a ). Thus, like the one-component case, 
the hopping term becomes 



i<j 

(19) 



and decoupled Hamiltonian is 



E 



(0) 



17 



!■/, -/ u ;S2T,92j ~ 2"(fiT +#U +ff2 T + 5 2l - 5) + 

+A 52 + 2(^11^21 +^11^21), (23) 



and the second-order perturbation energy is 

\{n\V\9)? 



M a > = EP , = V 

9 l9lT.ffl4.;92T.ff24.} 



P (0) p(0) 



E 



(2) 



(2) 



IT, 9 ^ ^U,9 ^ ^2T,g ^ ^21,g 



(2) 



(2) 



(24) 



(20) where (ct = -a) 
I 



E 



E. 



(2) 
i&,g 

(2) 

2cr,9 



.9io- + 1 



+ 



9la 



-Ugicj - + \± a U (g la - 1) + Fg^ - 

92a + 1 ff2<7 



-C/32 CT - Vg 2s + Ar - A U {g 2a - 1) + Vflto - ji a + A 



(25a) 
(25b) 
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Hence, the Landau expansion takes the form 
E 



The coefficients are 

b2a 



: -^{dli + 9n + 91] + 321 - 9) + Ag 2 + V (31T3U + 02f0u) - (MT3T + Ul9l) ■ 



92a 



1 



.92o- 



(26a) 
(26b) 

-Ug lr7 -Vg m ' U(g la - 1) + Vg lg - - p, a ' "' ^ 26 °^ 

To derive the stability condition, we need to consider the following derivative at the point Eq = 



-Ug 2 a-Vg 2 a+ fir- A 1/ (.92<t - 1) + Vg 2s 

gia + 1 gia 



C2<r 



A : 
= 2. 



0,V 2T = 0,^ 2i = 0): 










d 2 E n 


a 2 Bo 


a 2 Bo 


a 2 Bo 




d 2 E 


d 2 E 


a 2 Bo 


a 2 Bo 


D = 


di/ixidtpif 
d 2 E 


Wil 
d 2 E 


a 2 Bo 


a 2 Bo 




d 2 E 


d 2 E 


a 2 Bo 


dip2idip2i 

a 2 Bo 








dip2idtp2t 





2a 2 | 








2a 2 x 


C 2J 



C 2T 



26 2T 







2^21 



4a 2 i62j.c 2T - 4a 2 |62TC 2 j, + 16a2|a 2 j. &2T&2J.- 
I 



When V = 0, the insulating phase is no longer stable, 
thus a phase transition takes place. This critical condi- 
tion can be written as 4et2j&2T — c 2| or ^2^21 = c|^. 
(Either one being satisfied will destabilize the system.) 
As a result, there exists two phase boundaries, each one 



for a spin component, enclosing a region in the phase di- 
agram where one species is in the superfluid phase while 
the other one is still in the Mott insulating state. The 
phase boundaries for each can be written as (V = V/U) 



[(gi a -V) — p,„ + Vgw]{-gia ~\~ Mo - — ^9\s )[{jS2a - 1) - M<7 + A + Vg2a](-g2a + P*a - A - V g 2 &) 



(/V + 1 



V>Lcr)(/V 



Vg 2s ) 



(27) 



B. Phase diagrams 

Similar to the single component case, we first deter- 
mine the unperturbed groundstate and the range of the 
chemical potential, then depict the critical curve Eq. 
(|27|) . As V = V/U varies, there exists various filling con- 
figurations for the Mott state (unperturbed state). For 
the sake of clarity, we focus on the case V = 1, which is 
most accessible in current experiments. 

When the inter-species and intra-species repulsions 



are identical, the filling configuration for groundstate 
is very like the single component case, but there are 
much complicated spin-texture pattern, the spin imbal- 
ance m = g-\ — gi, (something like SDW) besides the 
CDW-like pattern in the former case. We list all the 
filling patterns for one spin component in the Appendix. 



After comparing E^ < E^i, E, 



AO) 



r(0) < E (0)^ we find 

the constraints on chemical potential \l a . For the most 
imbalanced case {g;0}, fi a € (g — l,g). For other filling 
patterns, the results are listed below. 



A 


g £ even 


g G odd 


ip- l,p) 


A*e (|[5 + (p-2)],|(ff-p) + A) 




(I[ ff -( P +1)] + A,l[.9+(P-1)]) 






£„€ (i[ 5 -(p + 2)]+A,i(.g+p)) 


/V G 


[l{g + ( p -l)],l[ g -(p+l)] + A] 





r 



With these, we plot the phase diagram for spin-| atoms in Figs. ([5]) and ([6]). The global structure of these dia- 
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0.4. 



0.3 - 



Superfluid 



Ml.{1:1} 



0.2 



Ml.{2:2} 
Ml,{3;3} 
Ml.{4:4} 
„, Ml.{5:5} 




2.5 



1.5 



0.5 




{4,0;0,3} 
{3,0;0,2|43,1;1,2} 



p.2;2,2} J 



3 5 



11,1,1,0}. ^ U3 . 30} (1,4,4,01 

r/ r i ^ 14.0:0,4} JfMf< 

i {3,0;0,3} ,p.l:l,3! 

,.....„ .{2,0;0,2| {±,i;i,i} ? ,{l,2;iJ f JfgVf 

!1,0;0,1} {1 l : l [J. {U;2,1J {1,3;3,1J Lj^f^ 

!0,1;1,0}I <0,2;2,0}> !0,3;3,0i r!",4;4,0i ^o',5;S,0} ' 

3 1 2 3 4 5 

fi/V 



10.3 



10.25 



H0.2 



0.15 



0.1 



0.05 



FIG. 5: (color online) The phase diagram spin-| atoms in two- 
component Bose-Hubbard model when V = U. MI stands for 
Mott insulator. gi±; g^\ , <72j.} reflects the filling config- 

uration in Mott state, where gi a particles with spin-cr are 
located in the Z-th site in a supercell. 



grams are similar to the one-component case, but there 
exist many small lobes hiding in a large lobe, as shown in 
Fig. (O . They are the phase boundaries of various fill- 
ing patterns in the Mott phase. The details are listed in 
Fig. where the sequence reflects the actual position 
of lobes in Figs. ([5]). For example, {1,0; 0,1} is writ- 
ten above {0,1; 1,0}, this means that in the 3D phase 
diagram, the lobe of filling {1, 0; 0, 1} is also located be- 
yond the lobe of {0, 1; 1, 0}. The illustrate these multiple 
lobes more clearly, in Fig. ([7]) we draw the cross-sections 
of the 3D phase diagram, A = (a) and A = U (b). The 
structures of other cross-sections are similar to these two. 
Every lobes from bottom to top corresponds to a filling 
configuration in Fig. © . From them, we find that there 
may exist weak superfluid (WSF) regions between two 
lobes that have the same particle number g. In WSF re- 
gions, one spin species becomes superfluid while the other 
type remains insulating. 



FIG. 6: (color online) The contour diagram for two- 
component case when V = U, revealing different CDW and 
SDW orders (spin imbalance) . 
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Ml 
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2 3 



FIG. 7: The cross-section for two-component case when 
V = U, showing the multiple lopes and possible weak su- 
perfluid regions (WSF). SF stands for superfluid. Detailed 
filling conditions can be found in Fig. (HJ). 



C. Spin dynamics in biased superlattices 

In addition to the phase transitions that we describe 
above, in the Mott insulating phase, the two-component 
BEC trapped biased superlattices can exhibit spin dy- 
namics and SDW-like patterns at different filling config- 
urations. 

At unit filling, Hubbard model [Eq. |T])] can be trans- 
formed into a Heisenberg-type model, via a canonical 
transformation and a second-order perturbation over fQ, 



namely, 

H cS = -jJ2(S-Sy + + ( J - J s ) J2 SIS*, (28) 

(«') <»i) 

where the (pseudo) spin operator is defined as usual, 

S? = 3 Yluu "lp cr A"^w with a a the Pauli matrix, a = 
x, y, z. The effective exchange energy is 



j ^ V j 

J ~ V 2 - A 2 ' s 



8t 2 U 
U 2 - A 2 ' 



(29) 
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When V < U and V > U, J = 0, then # cff = 
— JsYl/ij) ^iSji which is an Ising model. There is no 
spin wave excitations in this case. 

When V = U, J s = 2 J, then H cS = - J£ W ) S 4 • S,, 
which is a Heisenberg model. For ferromagnetic case 
(A < U, J > 0), we define Holstcin-Primakoff (HP) 
transformation as 

S + = VlS-tfcc, §~ = cWzS-tfc, S z = S-c f c, 

(30) 

where bosonic operator (c) creates (annihilates) spin 
deviations, 5 ± = S x ± iS' 3 ', £ is the total sp in. For lo w 
excited states, spin deviation is small, thus y/2S — &c « 
V2S; applying this transformation up to second order, 
the ferromagnetic Heisenberg Hamiltonian becomes 

H e s = —J 2^ &i ' 

<«) 

= E a + 2zJSj24 d i - JS J2(4 £ j + ft-c)(31) 

1 to') 

Here the ground state energy is E$ = —NzJS 2 with 
N the number of lattice site. This Hamiltonian is easy 
to diagnolize, the excitation energy in square lattice (in 
arbitrary dimensions) is 

e k = 2JS[z - cos (k • 6)} rj JS^ (k • 5) 2 = 2JSk 2 . 

S 8 

(32) 

Here, we set the lattice spacing to unity. 

For anti-ferromagnetic case A > U, J < 0, the 
single-occupied state is metastable and it would decay 
to a triplet groundstate which satisfies the Libe-Mattis 
theorem 1271 . However, this metastable state can be 
prepared [28j and probed 29] in an experiment-accessible 
time interval [l6j. if the band width of a single parti- 
cle t is much smaller than its band gap AE = A — U, 
t <C A_e[29|. Under this circumstance, we can perform 
the HP transformation in A-B sublattices and find the 
long-wave excitation in square lattices is 

e fe = 2z \J\ S^Jl^l t*y/2z\J\ Sk, (33) 

where jk = z^ 1 £ 5 cos (k • S). 



V. CONCLUSIONS 

In this paper, using a decoupling approximation, we 
analyze the possible phase diagrams of one-and two- 
component Bose-Hubbard models in optical superlattices 
in the mean-field level. As the potential bias A of the 
superlattice, the atomic repulsion U and hopping t, the 
filling configuration and chemical potential /j, varies, we 
discover complex phases in different parameter regions. 
For one-component case, there exists Mott states with 
CDW order and corresponding MI-SF transitions. For 
two-component case, besides the CDW in Mott state, 
there also exists weak SF regions where one spin com- 
ponent holds the CDW order while the other component 
becoming superfluid. In addition, the spin imbalance for 
a certain filling configuration (see Appendix) implies the 
existence of an SDW-like order in the Mott state. We also 
calculate spin dynamics of the two-component model at 
unit filling. The results explicit different low-energy dis- 
persions for different A/U. The features can be tested 
via many accessible probing techniques] 30, 31] in current 
experiments with single-site resolution [32|. 



Acknowledgments 

This research is supported by NCET, NFSC Grant No. 
10874017. 



APPENDIX A: FILLING CONFIGURATION FOR 
ONE-COMPONENT CASE 

The groundstate filling configurations {51552} for dif- 
ferent g (total number of particles in a supercell) and 
A = A/U are listed in the following table. Here, p is 
an even number, p — 1 > 0, reflecting the ratio between 
potential bias A and inter-atomic repulsion U. The no- 
tation {51; 52} means that there are g\ particles in the 
left (deep) site and 52 particles in the right (shallow) site 
in a supercell. 



A 


g £ even 


g £ odd 




r 9 1 p. 9 p\ 
1 2 ' 2 ' 2 2 / 

{5,0}, if g<p 


{£ji + (|-l);^-( 5 -l)} 

{ ff ,0}, if g<p-l 




{ B 1 P. 9 P\ 
\ 2 "1" 2 ' 2 2 J 

{g, 0} , if g < p 


f 9+1 1 p . 9-1 p \ 

X 2 2' 2 2 J 

{ ff ,0}, if g< P + l 



APPENDIX B: FILLING CONFIGURATION FOR p is still an even number that set the value of A, k,l,q 
TWO-COMPONENT CASE 



The filling configurations {gii, <?ij,; g2], 521} for two- 
component case is listed in the following tables. Here, 
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are integers that take values in certain intervals. The no- 
tation {gif, gii; <?2t) 92i} means that there are g^ spin-j 
atoms and g\i spin-J, atoms in the left (deep) site, gi-\ 



spin- 1 atoms and gi\ spin-J, atoms in the right (shallow) 
site in a supercell. 



A 


g G even 


(p- 


+ k - I - q,2 - {k - I) + q;l + q,Z - % - I - q} 

fee [o,f], le[o,k],qe fo,f-f-i] 


(p,P+l) 


{f + A: - Z - g, | - (fc - + g; / + q, § - § - I - q} 

fee [o,f],Ze[o,fc],ge [o,f-| -z] 



A 


(7 e odd 


(P- !,P) 


{a±l + (A ._ _ 9j (|_l)_ (fe _ + g; j + g ,a=l_(E_l)_ (/ + g) } 
fee [0,4i],/e[0,fc],ge [o,fi±i-£-i] 


(p,p+l) 


fce [o,Jti],/e[o,fc],ge [0,^-2-/1 



We can also dchnc the spin imbalance (mi = g^ — <?ij,m 2 = i?2T — <?2j.) for a given filling pattern, the results are 
listed in the following table. 



A 


g € even 


g S odd 


(p- 1,P) 


mi = § - f + 2 {k - I - q) , m 2 = 2 (/ + q) - (f - §) 


mi = 4^ - 1 + 2 (k - I - q) , m 2 = 2 (J + g) - - %) 


(p,P+l) 


mi = § - % + 2 (fc - I - q) , m 2 = 2 (/ + g) - (f - % ) 


mi = fi±i - | + 2 (fc - I - q) , m 2 = 2 (/ + g) - (2=1 - E) 



For filling configuration {3; 0}, we combine both filling pattern and spin imbalance in the following. 



A 


(/ e even 


<7 e odd 


(p- 1,P) 


{§ + fc,f-fc;0,0},.g<p,fce [0,f] 
mi = 2k, m 2 = 


{a±i + A5 ) fl=i-A!;0,0}, ff <p-l ) fcG [0,^] 
mi = 2k + l,rri2 = 


(P,P + 1) 


{f + fc,f-fc;0,0},.g<p,fce [0,f] 
mi = 2k, rri2 = 


{a±i + fc,a=i-A ! ;0,0}, fl <p+l ) fce [0,41] 

mi = 2fc + 1, m 2 = 
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